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Abstract 

We consider a mean-field dynamical urn model, defined by rules which give the 
rate at which a ball is drawn from an urn and put in another one, chosen amongst 
an assembly. At equilibrium, this model possesses a fluid and a condensed phase, 
separated by a critical line. We present an analytical study of the nonequilibrium 
properties of the fluctuating number of balls in a given urn, considering successively 
the temporal evolution of its distribution, of its two-time correlation and response 
functions, and of the associated fluctuation-dissipation ratio, both along the critical 
line and in the condensed phase. For well separated times the fluctuation-dissipa- 
tion ratio admits non-trivial limit values, both at criticality and in the condensed 
phase, which are universal quantities depending continuously on temperature. 
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1 Background: on dynamical urn models 



Dynamical urn models are simplified models of physical reality, which have always played 
an important role in the elucidation of conceptual problems of statistical mechanics and 
probability theory. 

The ancestor and prototype of this class of models is the Ehrenfest urn model [Q. It 
was devised by P. and T. Ehrenfest, in their attempt to critically review Boltzmann's 
if —theorem. Consider N balls, labeled from 1 to A^, which are distributed in two urns (or 
boxes). At each time step a ball is chosen at random (i.e., an integer between 1 and 
is chosen at random), and moved from the box in which it is, to the other box. Let 
A^i (respectively, A^2) be the numbers of balls in box number 1 (respectively, number 2) 
(A^i + A^2 = N). If the process is repeated indefinitely, for any initial condition the system 
will relax to equilibrium, characterized by a binomial distribution of balls in, say, box 
number 1: 

/.,eq = mi = /^)= (^)^- (1-1) 

This result is both intuitively translucent, and easy to derive (see below). The partition 
function of the system is equal to (for 2 boxes and A^ balls) 

nN 

Z{2,N) = —. (1.2) 

There are 2^ possibilities of distributing the balls amongst the 2 boxes, however all A^! 
labelings of the balls are equivalent. 

Finding the distribution of balls after n steps, fk{n) = V{Ni{n) = k), requires more 
effort. Kohlrausch and Schrodinger found the master equation for fkin), which they 
interpreted as the probability distribution of the position at time n of a random walker 
(played here by Ni{n)). The full solution of this master equation was given later on by 
Kac, Siegert and Hess 0, 1^, |§, §] . 

The Ehrenfest model is at the origin of a whole class of dynamical urn models, which 
we name the Ehrenfest class. They generalize the original Ehrenfest model in two ways, 
which we detail successively. 

The first generalization consists in considering M boxes instead of two. Then, at 
equilibrium, the joint distribution of the occupation numbers A"!, A^2, • • Nm, with 

M 

i=l 

is the multinomial distribution 
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as a simple reasoning shows. The marginal distribution of A^i is obtained by summation 
upon the other variables, and reads 

A,., = mi = *) = p)^. 

which is a simple generalization of (|1 . 1| ) . In the thermodynamic limit — oo, M — oo, 
with fixed density p = N/M, this yields a Poisson law: 

fi. = e~P— 

The partition function of the system is equal to (for M boxes and balls) 

Z{M,N) = ^, (1.3) 

as a simple generalization of the reasoning leading to (|1.2| ) shows. 

The dynamics of the Ehrenfest model and of its generalization to M boxes takes 
place at infinite temperature because there is no constraint on the move of the drawn 
ball. The second direction of generalization consists in defining these models at finite 
temperature, by introducing energy. We assume that the energy is a sum of contributions 
of independent boxes: 

M 

E(iVi,...,iVM) = E^(^0, 

i=l 

and choose a rule obeying detailed balance for the move of the drawn ball. For instance, 
for the Metropolis rule, the move is allowed with probability min(l, exp(— /^AE')). Heat- 
bath dynamics is another possible choice (see below). 

The fundamental ingredients for the definition of the models belonging to the Ehren- 
fest class are therefore 



(i) the statistics: a ball is chosen at random, and put in a box chosen at random, 

(ii) the choice of the energy of a box E{Ni), and of a dynamical rule (Metropolis or 

heat-bath), 

(iii) the geometry: for instance, boxes may be ordered on a line, or on the contrary be 
all connected. For short, we designate the latter geometry as the mean-field case. 

The backgammon model is a representative of the Ehrenfest class, corresponding 
to the choice (where 6 is Kronecker symbol) 

E{N,) = -6{N,,0). (1.4) 

This model has been extensively studied, mainly in its mean-field formulation p, ^ p!0| . 
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The dynamical (and equilibrium) properties of the Ehrenfest class depend crucially 
on the choice of statistics described in (i), which we will briefly refer to as the ball-box 
statistics. However, other choices are possible, which define new classes of dynamical 
urn models. The class which we will refer to for short as the Monkey class, because 
it corresponds to the image of a monkey playing at exchanging balls between boxes, is 
defined by: 

(i) the statistics: a box is chosen at random, from which any ball is drawn, and put in 

another box, chosen at random {box-box statistics), 

(ii) the choice of energy and dynamical rule (as above), 

(iii) the geometry (as above). 

A first example of a model belonging to this class corresponds to taking definition ( |1.4| ) 
for the energy |Q. This model, referred to as model B in 0, possesses non-trivial dy- 
namical properties |TU|. 



A second example, inspired from quantum gravity, corresponds to taking |TT[ 

E{N,) = \n{N, + 1). (1.5) 
It presents interesting properties both at equilibrium |Tl| and in nonequilibrium situ- 



ations [12]. In contrast with the backgammon model, or with model B, it possesses a 



phase transition between a fiuid phase and a condensed phase at finite temperature [|TT 
We will refer to this model as the zeta urn model, for reasons which will appear clear 
in the sequel. The present work is entirely devoted to the study of the nonequilibrium 
behavior of the zeta urn model in the mean-field geometry. 

Before specializing to this model, let us first present, in parallel, some formalism which 
applies to the two classes of models defined above, in order to underline the fundamental 
role played by the choice of statistics for both the equilibrium and nonequilibrium prop- 
erties of the models. Note that the equilibrium properties of the dynamical urn models 
defined above are independent of the geometry, because boxes are independent. 

For the Ehrenfest class, the partition function reads 

Z{M,N) = j:...j:^...^6hN,,N], (1.6) 

where 

is the unnormalized Boltzmann weight attached to box number i. For the Monkey class 
we have 

Z{M, N) = J2-..Y.p^,-.. p^,, 6 N„ n] . (1.7) 

Ni Nm \ i / 



4 



Using the integral representation 2i7r (5(m,n) = § dz z"^ " ^, we obtain 

dz 



where 



P{z) = E TT (Ehrenfest), (1.9) 

fc=o ^- 

oo 

= Y.Pk^'' (Monkey). (1.10) 

The equihbrium properties of the models are therefore entirely encoded in the tempera- 
ture-dependent generating series P{z). The presence or absence of the factorial term k\ 
has direct implication on the analytic structure of the series, and by consequence on the 
possible existence of a phase transition at finite temperature. 

The equilibrium probability distribution of the occupation number A^^i reads, for the 
Ehrenfest class, 

/fc,eq = V{N, = k) = {6{N,,k)) 

^E...E^(-..)||...^.(E-.iv) 



Pk Z{M-l,N-k) 
k\ Z{M,N) 

For the Monkey class one obtains 



:i.ii) 



Z{M -1,N - k) 

At infinite temperature, equation (p..8|), together with ( |1.9| ) and (|1.10|) , respectively lead 
to ( |1.3| ) and to 

Z(M.iV)^'f^ + ^-^>'. (1.13) 

^ ' ^ (M-l)!iV! ^ ' 

In the thermodynamic limit at fixed density p, the free energy per box is defined as 

pF = - Urn ^lnZ{M,N), N ^ Mp. 

At infinite temperature, equations ( |1.3| ) and ( |1.13| ) yield 

lim 13 F = plnp — p (Ehrenfest), 

= plnp-(p + l)ln(p + l) (Monkey). 

At finite temperature, the free energy can be obtained by evaluating the contour integral 
in ( |1.8|) by the saddle-point method. The saddle-point value Zg of 2; is a function of 
temperature and density through the equation 
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and the free energy per box reads 

pF = p\nZs-\nP{zs). (1.15) 

Similarly, we obtain the following expressions for the equilibrium occupation probabilities 
in the thermodynamic limit 

k 

/Mq = Ifp^ (Ehrenfest), (1.16) 

k 

= Pkj^ (Monkey). (1.17) 

This formalism will be illustrated below on the zeta urn model, studied in this work. 
For both classes of models, the temporal evolution of the occupation probabihty 

f,{t)=V{N,{t) = k) 

is given by the master equation 

dfkit) 



, ^{^^k+i/ + 'n'£,k-i - T^k,e - 'n'e,k) , (1-18) 

where vr^, ^ denotes the contribution of a move from the departure box, containing k balls, 
to the arrival box, containing i balls. Restricting our study to the mean-field case, we 
have (for k,i > 0) 

T^k/ = kfkfeWk/il - 4,o) (Ehrenfest), 

= fkfeWkA^-6,,o) (Monkey), (1.19) 

where the term 1 — 6k,o accounts for the fact that the departure box cannot be empty, and 
where the acceptance rate Wk/ depends on the dynamics chosen. With the Metropolis 
rule we have 

Pk-m\ 



Wkj = min 1, 



PkPe+i I 



while the heat-bath rule leads to [T^ 



W^M=^ E/.^ , (1-20) 
Pi V£=o Pi J 

which only depends on the label i of the arrival box. 

In all cases, equation ( |1.1^ ) can be seen as the master equation of a random walk for 
A^^i, i.e., over the positive integers A; = 0, 1, . . ., 

^•^''^^^ = fik+i fk+i + h-i fk-i - {f^k + h) fk (fc > 1), (1-21) 



dt 

dm 
dt 



/^i/i-Ao/o, (1.22) 
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generalizing the result of Kohlrauscli and Schrodinger for the Ehrenfest model. In these 
equations, 

M - f^k - 2^ —r- 

1=0 Jk e=o Jk 

are, respectively, the hopping rate to the right, corresponding to A^i = A; ^ A^^i = A; + 1, 
and to the left, corresponding to Ni = k ^ Ni = k — 1. The equation for /o is special 
because one cannot select an empty box as a departure box, nor can Ni be negative, 
hence A_i = fio = 0. In other words a partially absorbing barrier is present at site k = 0. 
The random walk is locally biased, to the right or to the left, according to whether its 
velocity Afc — fik is positive or negative, respectively. It is easy to verify that ( |1.21|) 
and ( |1.22| ) preserve the sum rules 

E^W = 1' (1-23) 

k 

Y.kfk(t) = {N,{t))=p, (1.24) 

k 

expressing respectively the conservation of probability and of the number of particles. 

The equilibrium occupation probabilities (|1.16| ), (|1.171) are recovered as the unique 
stationary state (d/^/dt = 0) of the master equations ( p..21j ), ( |1.22[ ). In contrast, except 
at infinite temperature, where Wk/ = 1, so that the rates and fik simplify, the master 
equations cannot be solved explicitly. The difficulty comes from the fact that the rates 
Afc and Hk are functions of the fk, hence the master equations are non-linear. However, 
as will be illustrated by the present work, the long-time behavior of these equations is 
amenable to analytic computations. 

To summarize, for the two classes of dynamical urn models described above, finding 
the properties of their equilibrium states is in general easy, but attaining dynamical 
properties is much more difficult to achieve. However if one restricts the study to the 
fluctuating number of balls in a given box, denoted by Ni{t) all throughout this paper, it 
is, in some cases, possible to predict the long-time behavior of its probabihty distribution, 
fk (t) , and also of its two-time correlation and response functions. This has been done in 
previous studies for the backgammon model P|, pUj]. 

In the present work we focus our interest on the mean-field dynamical urn model 
defined by the choice of energy (|1.5| ), and box-box statistics -the zeta urn model. We 
pursue the investigation of its nonequilibrium properties, initiated in ||l2l. The next 



section is devoted to a more complete presentation of the model and to an outline of this 
paper. 
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2 The zeta urn model 



At equilibrium, the zeta urn model is defined by its partition function ( |1.7| ), where the 
Boltzmann weight 

= e-^^W) = (iV, + 1)-^ (2.1) 

corresponds to the choice of energy ( |1.5|) ||11|| . The equilibrium phase diagram of the 
model |lT[ easily follows from the analysis of the previous section (see ( |1.8| )-( prT^ )). 
At low enough temperature (/? > 2), there is a finite critical density: 

_ p'ji) _ ap-i)-m 
p(i) m 

where C, denotes Riemann's zeta function. 

In the fluid phase (p < pc), the equilibrium distribution ( |1.17| ) decays exponentially, 
since < 1. 

At the critical density (p = pc), we have Zs = I, hence 

f - ^ - (k + l)-" 

which is known as the zeta distribution. In the regular part of the critical line (/5 > 3), 
the mean squared population is finite, and equal to 



k=0 



C(/?) 



while it is infinite in the strong-fluctuation case {2 < (3 < 3). Throughout the following, 
we shall restrict the study to the regular part of the critical line. 

In the condensed phase (p > pc), a macroscopic condensate of particles appears. 
Indeed, equation ( |2.2| ) still applies to all the boxes but one, in which an extensive number 
of particles, of order — Mpc = M{p — pc), is condensed. 

The dynamical definition of the model was given in the previous section. For heat- 
bath dynamics, equation ( |1.18|) , together with ( |1.19| ) and (|1.20|) , leads to [0 

dm 



dt 

dm 
dt 



with 



and 



/fc+i(t) + ait)rk-ifkm - (1 + (^it)rk) m (k > 1), 
/i(t)-a(t)ro/o(t), (2.3) 



Tk- - - ' ' 



Pk fk,eq \^ + 2, 

1 - m 



a{t) 



k=0 
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We assume that at time t = the system is quenched from its infinite-temperature 
equihbrium state to a finite temperature T = Hence, by ( p..l4| ) and ( |1.17] ), the 

initial occupation probabihties read 

In the fiuid phase (p < pc), the equihbrium distribution /fc^eq ( |1 ■ 1 7|) is a stationary 
solution of ( p.3|) , corresponding to cXeq = Zs- The convergence of fk(t) towards fk,eq is 
characterized by a finite relaxation time, depending on j3 and p. 

The long-time behavior of the distribution fk{t), both at criticality (p = pc) and in 
the condensed phase (p > pc), is the subject of the next section. In section 4 we establish 
the dynamical equations obeyed by the two-time correlation and response functions of 
the population of a given box. It is indeed well-known that nonequilibrium properties are 



more fully revealed by two-time observables [1^, |14[. Section 5 is devoted to the analysis 



of the equilibrium properties of these functions at criticality. Section 6 is devoted to the 
analysis of their nonequilibrium properties at criticality, including the violation of the 
fiuctuation-dissipation theorem. In section 7 we briefiy investigate the nonequilibrium 
properties of the model in the condensed phase. 



3 Long-time behavior of occupation probabilities 

3.1 At criticality {p = pc) 

We investigate how, starting from the disordered initial condition ( p.4|) , the occupation 
probabilities fkit) converge toward their critical equilibrium values /fc,eq, given by ( |2.2|) , 
which are the stationary solutions of equations ( |2.3| ) corresponding to a^q = 1. In analogy 
with the analysis done in [^], we anticipate that a{t) converges to this value as a power 
law: 

and consider two regimes: 

Regime I: k fixed and t ^ 1 

This is the "short-distance" regime, considering k as the position of a fictitious random 
walker, as in section 1. It is therefore analogous to the Porod regime for phase-ordering 
systems [|T^, . 

Setting 

/fc(t)~/fc,eq(l+^fet-"), (3.1) 

equation ( |2.3| ) yields 

Vk = Vo + Ak, (3.2) 
where vq and A are determined below. 
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Regime II: k and t simultaneously large (scaling regime) 

In this regime, we look for a similarity solution to (p. 3D of the form 



fk{t)^fk,e^F{u), u = kt-'/^. (3.3) 



The structure of the master equations ( |2.3| ) indeed dictates that the scaling variable is 
the combination kt'^^"^ . Starting from a random initial condition, for a large but finite 
time t, and for k much smaller than an ordering size of order t^/^, the system looks critical, 
i.e., the distribution /^(t) has essentially converged toward the equilibrium distribution 
/fc^eq- This implies -F(O) = 1. To the contrary, for k ^ t^/^, the system still looks 
disordered, i.e., the fkif) fall off very fast. Hence F(oo) = 0. It will indeed be shown 
below that fkif) ~ exp(— /c^/(4t)), as the scaling function falls of very fast for u ^ 1: 
F{u) ~ exp(— ^^74). Note the close analogy between the present situation and critical 
coarsening for ferromagnetic spin systems []T^ |r7|, where the scaling variable is |r| t^^/^, 
and where the role of fk{t) is played by the equal-time correlation function C(|r| ,t). 

In order to determine the exponent tu, we use the sum rules (|1.23| ) and ( |1.24|) , which 
yield respectively 

t--{vo + Ap,) = /i = — ^ / u-^{l- F{u)) d«, (3.4) 

C(/3) Jo 

t-^iVoPc + A/Xe) = r(^-^)/'/2, h = -7^J U'~^ (1 - F{u)) dw. (3.5) 



These equations are compatible only if the right-hand side of equation ( |3.4| ) is subleading. 
In the case of a regular critical point (/3 > 3), we thus obtain 

/?-2 1 

u) = > - 

2 2 



and 



A = ^, vo 



Pc- pV Pc- pT 



Inserting the form (p.3|) into ( |2.3|) , and using the fact that u > 1/2, leads to the 
differential equation 

VF{u) = 0, (3.6) 
where V is the linear differential operator 



The solution of (|3.6| ) is 

2-/3 



2 



Fin) = ^;^l y^e-y/'dy. (3.8) 
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We present in Appendix A an alternative way of solving equation (3^), using the 
Mellin transformation. In particular ( |A.2| ) yields the explicit expression 

Mi_F{f3 - 2) _ TT^^ 2''" 

C(/3) "(/3-2)r(^)C(/?)' 



As an illustration of the above, let us determine how the variance of the population 
of box number 1, Var A^i(t) = (A^i(t)^) - {Ni{t)f = {Ni{t)^) - pi, converges at long 
times to its equilibrium value pc~ pI- Only the scaling regime matters for the long-time 
behavior of this quantity, and of all the quantities at criticality to be considered hereafter. 
Using (|3.3| ), we obtain 

VarA^l(^)-(/.,-p^)^-A^-(^-3)/^ 
Mi_f(/3 - 3) 23-/3 



with, by ( ITOD , 



A 



CW) (/?-3)r^ c(/5) 



3.2 In the condensed phase {p > pc) 

We set 



and consider the same two regimes as at criticality (see |12] for more details). 

Regime I: k fixed and t ^ 1 (short-distance regime) 

Equations ( |3.1| ) and ( p.2| ) still hold, but now with = 1/2, for any [3 > 2. 

Regime II: k and t simultaneously large (scaling regime) 

Looking for a similarity solution of equations ( p.3| ), of the form 

/.(t)^^, u = kt-'l\ (3.9) 
we obtain for the scaling function F{u) the linear differential equation 



Following [|l^, we notice that the amplitude A is determined by the fact that equa- 
tion (|3.1CI|) has an acceptable solution F{u), vanishing as m ^ and u — > oo. The 
normalization of the solution F{u) is determined by the sum rule ( |1.24|) , which yields 



u F{u) du = p — Pc- (3-11) 

The parameter vq entering Regime I is determined by the sum rule ( 1.23| ), leading to 

Vo = -Apc- / F{u) du. 
Jo 
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At variance with equation ( p.6|) , the differential equation (|3.10 ) cannot be solved in 
closed form. It can be recast in Schwarzian form, without first derivative, by setting 



F{u) = uY{u)H{u), 

with 

Y{u)=u-^/^-^e^^/^-^"/'. (3.12) 
We thus obtain for H{u) a differential equation of the form 

nH = 0, (3.13) 

with ^ 

n = -^ + Wiu), (3.14) 

and where the potential W{u) reads 

Equation ( p.l3| ) is a biconfluent Heun equation [|1^ . We will therefore refer to H as the 
Heun operator, denoting its discrete eigenvalues by En and the corresponding eigenfunc- 
tions bu Hn{u). Equation (|3.13|) implies that the ground-state eigenvalue reads Eq = 0, 
while the associated eigenfunction Hq{u) is simply related to the scaling function F{u): 

F{u) = cuY{u)Hoiu). (3.16) 

By ( |3.11| ) we have 

. (3.17) 



oo 

u'^Y(u) Hn(u) du 



The spectrum of the operator H, and related quantities such as the scaling function 
F{u), can be further investigated in the limiting regimes of high and low temperature. 

High temperature {f3 0) 

The analysis of this limiting situation will be helpful in section 7, although it is 
of no direct physical relevance, since the condensed phase only exists at low enough 
temperature (/? > 2). 

For /5 = and A = 0, the potential (|315|) becomes W{u) = - 3/4, so that H 

is the Hamiltonian of a harmonic oscillator on the half-line m > 0, up to a scale, with 
Dirichlet boundary condition at u = 0. Its spectrum is En = n {n = 0,1, .. .), and the 
(unnormalized) ground-state eigenfunction reads 



i7o(M) =Me-" /^ (3.18) 

hence 

F{u) = ^ue--'/\ (3.19) 
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Perturbation theory can then be used to determine the small-/? behavior of the am- 
plitude A. Expressing that the lowest eigenvalue reads Eq = with no correction, we 
obtain 

2P' 



^-Au + Hl{u) dw + ■ ■ ■ = 0, 

hence 

A ^ -—j3. (3.20) 



Low temperature [13 oo) 

In this other limiting situation, the operator Ti simplifies as follows. If we rescale A 
and u according to A ~ m ~ all terms in the expression (|3.15| ) for the potential 

W{u) scale as (3. In other words, (3 plays the role of 1/h in Quantum Mechanics, and 
/? — ^ oo is the semi-classical regime. Expressing that the minimum of the potential is at 
zero yields the first estimates A'^ u (2/?)^/^. 

A more refined analysis consists in expanding the potential W{u) around its minimum. 
Setting 

A={2(3f/^-a, u = {2(3f/^ + 13^/^, (3.21) 
the operator Ti becomes 

«~+r'"(-il + 50,„'). (3.22) 

The expression in the parentheses is again proportional to the Hamiltonian of a harmonic 
oscillator. Expressing that the lowest eigenvalue of the right-hand-side of ( p.22 ) is Eq = Q 
determines a = 2^/^ — 2~^^'^P~^^'^, hence 

A ^(2/3)1/2 -2^2 + 2-1/2^3-1/4. 

The spectrum of Ti. reads En ~ n {n = 0,1, . . .), and the (unnormalized) ground- 
state eigenfunction reads Hq{u) ~ e"""^/"^, i.e., 

exp(-(^z|§!^!)!). (3.23) 

Both this expression and the function Y{u) become singular in the /? — >■ oo limit, so that 
a direct analysis of equation ( p.lO| ) is needed in order to derive the behavior of the scaling 
function F{u) at low temperature. The above analysis suggests to look for a solution 
depending on the scaling variable y = {2(3)'^^'^u. The term involving the second-order 
derivative d^F/dv? is then negligible. The simplified form of equation (|3.1CI|) , namely 

d P 

y{y -lf^ + {2y' - 1)F = 0, y= m-''^u, 
dy 
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admits the normalized solution 



where Ei is the first exponential integral {Ei{l) = 0.219383934). The scaling function 
F(u) is therefore nonzero in the /? — >■ cx) limit only for y < 1, i.e., u < (2/5)^/^. The upper 
bound y = 1, i.e., u = (2/5)^/^, coincides with the point where the eigenfunction Ho{u), 
as given by equation ( p.23| ), peaks at low temperature. 



4 Two-time observables: dynamical equations 

In this section we will successively establish dynamical equations for the two-time correla- 
tion function C(t, s), for its derivative dC(t, s)/ds, and for the response function R(t, s). 
(See for similar techniques.) 

We consider the (connected) two-time correlation function C (t, s) between the popu- 
lation of box number 1 at times s (waiting time) and t (observation time), with < s < t: 

C{t,s) = {Nr{t)N^{s)) - {N,{t)) {N,{s)) = {N,{t)N,{s)) - p'. 

This definition can be recast as [|l^ 

oo 
k=l 



with 



7fc(t, s) = Y.3 fA-^) HNiit) = k\N,{s)= j}. 



The evolution of the 7fc(t, s) with respect to t is given, for t > s, by a master equation 
similar to (121 



d'y (t ts) 

— ^ — = 7fc+i(t, s) + (T(t)rfc_i7fe_i(t, s) - (1 + cr{t)rk) 7fc(t, s) {k>l), 
a7o(t,s 



dt 



7i(t,s)-a(t)ro7o(t,s). (4.1) 



These equations preserve the sum rule J2klk(t,s) = p. At t = s, the initial conditions 
are 7^(5, s) = kfk{s), implying 

00 

C{s,s) = Y.k'fk{s)-p\ 

k=l 

which is the variance of the population of box number 1 at time s, as it should. 
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In the discussion of the fluctuation-dissipation theorem, we will need expressions of 
the time derivative dC(t, s)/ds. We have 

— Q = 2^kipk{t,s), 

k=l 

with 

<97fc(t,s) 



ds 

The evolution of the fk(t, s) with respect to t is again given, for t > s, by equations 
similar to (|2l 



— ^ — = fk+i{t, s) + a{t)rk-i^k-i{t, s) - (1 + (T(t)rfc) Lpk{t, s) (A; > 1), 
ot 



ipi{t,s) - a{t)ro(po{t,s), (4.2) 



dt 

with initial conditions at t = s 

iPk{s, s) = -/fc+i(s) + a(s)rfc_i/A:_i(s) {k > 1), 

so that 

The two-time response function R{t, s) is a measure of the change in the mean popu- 
lation of box number 1 at time t, induced by an infinitesimal modulation of the conjugate 
variable, i.e., the local chemical potential acting on the same box, at the earlier time s. 

In the presence of an arbitrary local, time-dependent chemical potential the 
energy of box number 1 at time t reads 

E (iVi(t)) = ln(iVi(t) + 1) - ^(t)iVi(t). 

The occupation probabilities of this box now depend on fi{t): we denote them by fk{t). In 
the thermodynamic limit, i.e., to leading order as M — ^ 00, the occupation probabilities 
of all the other boxes (z = 2, . . . , M) are still given by the fk{t)- 
The response function reads 



IX=Q k=l 

with 

hk{t,s) 



fi=0 
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The modified occupation probabilities fj^it) obey the dynamical equations 

= + ^(t) e^^(*V,_i/^,(t) - (l + a{t) e/^'^Wr,) /^t) > 1), 

= - e^-Wro/o^t). (4.3) 

The initial values (0) = /fc(0) and the parameter a{t) are unchanged. The dynamical 
equations ( [4.3| ) preserve the sum rule /fc (^) = 1- 

Equations ( |4.3| ) imply that the hk(t, s) obey, for t > s, 

dhk{t, s) ^ ^^^^(^^^ _j_ (x{t)rk-ihk-i{t, s) - (1 + a(t)rfe) /;,fc(t, s) (/c > 1), 



dt 
dhoit, s 



hit,s)-ait)rohoit,s), (4.4) 



dt 

with initial conditions at t = s 

hk{s, s) = (3a{s) (rfc_i/fe_i(s) - rkfk{s)) {k > 1), 
ho{s,s) = -/?a(s)ro/o(s), 

so that 

R{s,s)=/3{l-fois)). 

The behavior of the two-time observables will now be successively investigated in 
the next three sections, first at criticality (p = pc), both at equilibrium and in the 
nonequilibrium regime, and then in the condensed phase (p > pc). 



5 Equilibrium critical dynamics 

When the waiting time s becomes very large, keeping the difference r = t — s > fixed, 
two-time quantities reach their equilibrium values, which only depend on r, both in the 
fluid phase (p < pc) and along the critical line (p = pc). This section is devoted to the 
latter case. 

The equilibrium correlation function Ceqir) reads 



Ceq(r) = ^/i;7fc,cq(r) - pi 



k=l 



where the equilibrium values 7A:,eq(T) of the 7fc(t, s) obey 

d7fc,eq(^) 



dt 

d7o,eq(r) 

dt 



7fc+l,eq(r) + rfc_i7fc_i,eq(r) - {1 + Tk) 7fc,eq(^) (^ > 1), 

7l,eq(r) -ro7o,eq(T), (5.1) 
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with initial conditions 7A;,eq(0) = fc/fc,eq- This imphes consistently that Ceq(O) = jic — p1 
is the equilibrium variance of the population of a generic box. 

As the t ime difference t gets large, we have 7fc,eq(''") — pc fk,eqi so that Ccq(T) — 0. 
The decay of Ceqir) for large r can be investigated along the lines of section 3.1. The 
two regimes I and II are again to be considered separately, although results concerning 
the former will not be needed explicitly. In the scaling regime, we look for a similarity 
solution to equations ( ^.1|) of the form 

7fc,cq(r) ^/fe.eqrl/^GeqH, U = kT-^/^ (5.2) 

for which equations ( |5.1| ) yield the differential equation 

'^ + 1) Gcqi^) = 0, (5.3) 



2. 

with boundary conditions Geq(O) = and Geq{u) ^ u as u oo. The differential 
operator V was defined in equation (|3.71). 

Equation ( p.3| ) can be solved by the method of Appendix A. With the definition ( |A.1| ), 
we have the functional equation 

MgJz + 2) _ z-1 
MgJz) 2{z + 2){P-1-z)' 

whose suitably normalized solution reads 

The decay of the equilibrium correlation function at large r is again dominated by 
the scaling regime, and the scaling form ( |5.2|) yields 

Ceq(r)^Aeqr-(^«-3)/2, (5.4) 

with 

MGjf3-2) vrr(/^) 

As expected, the fluctuation-dissipation theorem holds at equilibrium. Indeed the 
equilibrium values (pk,eq{'r) and /ifc,eq(T) of ipk{t, s) and hk(t, s) obey the same equations, 
identical to equations ( |5.1| ), with initial conditions: 

^fc,eq(0) = P(pk,eq{0) = (3 {fk,eq " /fc+l,eq) (^ > 1), 

/iO,eq(0) = /5v50,eq(0) = -/?/l,eq, 

hence the identities 

hk,eq{r) = P'fk,eq{r) (5.5) 

and 

The last formula is the fluctuation-dissipation theorem in its usual differential form. 
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6 Nonequilibrium critical dynamics 

We now turn to the nonequilibrium critical behavior of the two-time correlation function 
C(t,s), response function R(t,s), and fluctuation-dissipation ratio 

R{t,s) 

^(M) = ^^^, (6.1) 
ds 

in the scaling regime where both time variables s and t are large and comparable. Here- 
after X will denote the dimensionless time ratio 

t 

X = - > 1. 
s 

We first analyze the scaling behavior of the correlation function C{t, s). Looking for 
a two-variable scaling solution to equations ( [4.1|) of the form 

-fk{t,s)^ fk,ecit'/^G{u,x), u = kt-'/\ x = t/s, (6.2) 

we obtain the partial differential equation 

(^x^ + V + ^'jG{u,x) = (6.3) 

for the scaling function G{u,x), with initial condition 

G{u,l) =uF{u), (6.4) 

and boundary condition G{oo,x) = for all x > 1. The function F{u) is known 
from (|3.8|), and the operator V is given by (|3.7]). 

As shown in Appendix B, equation ( |6.3D can be solved explicitly by the method of 
separation of variables. Equation (|6^ ) imphes the scaling law 

G{t,s)^s-^^-^^/^^x), (6.5) 

where 

m Jo 



Hx)= u'-^G{u,x)du. 



Using (|B.4|) , we obtain 

$(a;) =x-^/2^A„x-", (6.6) 



oo 



n=0 



with 

An 



a„MG„(/j-2) 

m 
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More explicitly, equations ( B.3 ) and ( [B.5|) imply that the leading coefficient of the ex- 
pansion ( |6.6D reads 

^1/2 23-/3rfm) 

= , (6.7) 

3(/5 + i)r(^) c(/3) 

while the coefficient ratios are rational functions of (3: 

Ao ~ 2- -LV^^ + ■ (2A; + 3)fc!(n - A;)! /J, 2£ + /3 + 3 ' ^^'^^ 



I.e., 

/3(2/3 + 3) A2 /3(/3 + 2)(8/32 + 52/3 + 45) 



etc. 



Ao 10(/? + 3)' Ao 280(/5 + 3)(/?+5) 

The expansion ( |6.6| ) is convergent over the whole physical domain (x > 1). For x — > 1, 



i.e., r = t — s <C s, the equilibrium result ( |5.4| ) is recovered as 

$(x) ^ v4eq(a; - l)-(/3-3)/2_ 

These properties can be checked by noticing that the expression (|B.lj ) for the Laguerre 
polynomials simplifies to 

/9„l/2\ (/3+l)/2 

4^)72(^2/4) ^ j^^j J(,^.i)/2(nnV2), 

where J is the Bessel function, in the scaling regime where the order n is large and u is 
small. The subsequent integrals can be estimated for n large, yielding 

A ^ ^ „(/3-5)/2 



2^r(f) c(/5) 



This asymptotic expression ensures the convergence of the series (|6.6| ) for all a; > 1, and 
estabhshes (p.9[), including the prefactor. 

We now turn to the derivative dC{t, s)/ds of the correlation function. Looking for a 
two-variable scaling solution to equations ( [4.2| ) in the scaling regime, of the form 

Mt,s)^ fk,e^r'/^G^'\u,x), (6.10) 

we get the partial differential equation 

(^x^ + V-^^G^'\u,x) = 0, (6.11) 

with initial condition 

G«(m,1) =/3^^-2F'(m). (6.12) 
u 
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Equation ( 5.11 ) can again be solved by the method of Appendix B. We thus obtain 



ds 

where 



-(/3-l)/2 

= — / u^'^G'^^\u,x)du. (6.14) 



C(/3) Jo 

This expression will not be made more explicit, as equation ( |6.5|) yields more directly 

^WU) = ^^^x) ~x^, (6.15) 
2 ax 

i.e., 

oo 

$(i)(a;) = 5] A« (6.16) 

n=0 

with 

4'^ = (6.17) 

The scaling behavior of the response function R{t, s) can be determined by the same 
approach. Looking for a two-variable scaling solution to equations ([4.4|) in the scaling 
regime, of the form 



equations (f4.4|) yield the partial differential equation 

(^x-^ + V-^^G^'\u,x) = 0, (6.18) 

with initial condition 

G^^\u,l) = p''^-/3F'iu). (6.19) 
u 

We thus obtain 

R{t,s)^s~^^~^^/^^^^\x), (6.20) 

where 

<l>(^)(x) = ^——— u^~^G^^\u,x)du, (6.21) 



a/3) Jo 
i.e., explicitly 

$(2)(a;)=a;-/3/2^42)^-n_ (g^22) 
n=0 

The leading coefficient of this expansion reads 
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while the coefficient ratios are again rational functions of (3\ 

A^^^ 1 n rg ■ I " {-f{2k + p + i) ^^ 2^ + p + 2 

i.e., 

^ ^ (3{2(3 + 1) ^ /j(/3 + 2)(8/j^ + 28/3 + 9) 

42)-6(/?+l)' 42)- i20(/3+l)(/3 + 3) ' 

The scaling results and ( |6.2(J| ) imply that the fluctuation-dissipation ratio 

X(t, s), defined in equation (|6.1|) , only depends on the time ratio x in the nonequilibrium 
scaling regime. We have indeed 

The scaling function X[x) turns out to be universal, whereas $*^^)(a;), and 

which respectively enter equations (|b.5|), ( |6.13| ), ( |b.2(JD , are only universal up to a scale 
fixing. (A definition of universal quantities in the present context will be recalled in the 
Discussion.) 

Equations (|6.16| ) and (|6.22|) yield the expression 

00 

X{x) = = E (6.24) 

n=0 

which is again convergent for a; > 1. The limit fluctuation-dissipation ratio, X^o = 
X{oo) = C,o, takes the simple value 

X. = >1±1. (6.25) 

Equation ( |6.24D also yields 

_ /j^(-2/3=^ + 2/3^ + 78/3 +117) 
3(/3 + 2)(/? + 3)' 45(/? + 2)(/5 + 3)2(/5 + 5) ' 

The behavior of the fluctuation-dissipation ratio X{x) close to equilibrium (i.e., for 
X — > 1) is also of interest. Equations (|6.14| ) and (|6.21|) for a: = 1, together with ( |6.12|) 
and ( |6.19| ), imply 

lim ((3^^^\x) - <l>(2)(x)) = --L^MF,{f3-2) '^^ ^ ^ 



This expression, together with equations (|6.9|) and (|6.15|) , yields 

2rf^f 

X{x\ ^ 1 ^ (x - l)(^-i)/2 
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confirming that the fluctuation-dissipation theorem is restored for x — 1, and exphcitly 
giving the leading violation of that theorem in the scaling regime, which is proportional 
to (r/s)(^-i)/2. 



7 Nonequilibrium dynamics in the condensed phase 

We finally turn to the long-time behavior of the two-time quantities C{t, s), R{t, s), and 
X(t,s) in the condensed phase (p > pc). 



In Regime II, we look for two-variable scaling solutions to equations (|4.1| ) 
and (|4.4| ), inspired by (|3l9|) , of the form 



7fc(t,s)^^G(n,x), Mt,s)^^G^'\u,x), h,{t, s) ^ ^G^'\u,x), (7.1) 
t 

where u = kt~^/'^ and x = t/s. 

The three scaling functions obey partial differential equations of the form ( |6.3|) , ( |6.11| ), 
or ( |6.18| ), where the differential operator V is replaced by 



^ f u ^ /3 + 2\ d 3 A 



This differential operator is related to the Heun operator Ti. ( ^.141 ) and to the function 
Y{u) defined in equation (|3.12| ) by the conjugation 



Y{u)n = VY{u). (7.3) 

Hence V and 7i share the same eigenvalues En {n = 0,1,...), which are not known 
explicitly, as already mentioned, except Eq = 0. 

The two-time scaling functions defined in ( [7.1|) are then determined in analogy with 
the critical case. We thus obtain 

OD OO 

G{u,x) = Y{u)Y.anX-''--^'^ Hniu), G^^^^\u,x) = F(m) ^ ^i'^)^-^"^^/' //„(«), 

n=0 n=0 

(7.4) 

where the coefficients determined by the initial conditions 



G{u,l)=F{u), G^'\u,1)=(--^]f{ 



u)--F'(u), G^^UuA) = -^F'(u). 



u u 

(7.5) 

Using the expansions ( [7.4|) , as well as the linear dependence of F{u) in (p — p^) (see 
equations ( p.l6|) , ( p.l7D ), we obtain the following scaling forms for the correlation and 
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response functions: 

oo 

Cit, S) ^ S'/^ (p - pe) HX), $(X) = E 

n=0 

oo 

-'/^p-p,)<^^'\x), = E^'^^"""", (7.6) 



s 



R{t, S) ^ S-^/^ (p - pe) <l>(2)(x), = E 4'^ 



n=0 
oo 



n=0 



The scaling functions only depend on /?. They take finite limit values, 

both at X = and at x = oo. The fluctuation-dissipation ratio 

X{t,s)^X{x) ^ ' 



/3<l>«(x) 



again only depends on the time ratio x in the nonequilibrium scaling regime. The whole 
scaling function X{x) is universal, just as in the critical case. In particular, the limit fluc- 
tuation-dissipation ratio Xoo = X{oo) = A'q^ /Aq'^ in the condensed phase is universal, 
and only depends on /?. Using equations ( |3.16| ) and ([7.5|), we can derive the following 
expression for X^o in terms of the ground-state eigenfunction Hq{u): 

^^^h^^^^^ . (7.7) 

u Hq{u) du 



The above expressions cannot be made more explicit in general. More quantitative 
results can be derived at high temperature and at low temperature, using results of 
section 3.2. 

Let us consider the limit fluctuation-dissipation ratio Xoo, as given by equation (|7.7|). 
In the high-temperature case, by (|3.18|) and ( p.20|) , we obtain 

X^ = ^-in-2)P + --- (7.8) 

In the low-temperature case, using ( p.21|) and (|3.23| ), we obtain 

Xoo = r'/'-^ + --- (/5-00). (7.9) 

It is instructive to compare the limit fluctuation-dissipation ratios corresponding to 
the critical point and to the condensed phase. The value of X^o at criticality is given by 
the analytical expression ( |6.25| ). The value of Xoo in the condensed phase is obtained 
by a numerical evaluation of ( |7.7| ): the ground-state eigenfunction Hq{u) is obtained by 
numerically solving the differential equation ( |3.10| ) or ( p.l3| ), and then used to evaluate 
the integrals entering the result (|7.7|) . The values thus obtained smoothly interpolate 
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between the limiting laws (|7.8| ) and In Figure |T], thick full lines show the physical 

values of the fluctuation-dissipation ratios, while thin dashed lines show their continua- 
tion to high temperatures. Both fluctuation-dissipation ratios start from 1/2 at infinite 
temperature, and converge to the limit values Xqo = and X^o = 1 at zero temperature, 
respectively in the condensed phase and at criticality (see Discussion). 

1.0 



0.8 



0.6 



0.4 



0.2 




4 
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Fi gure li Plot of the limit fluctuation-dissipation ratio X^q against inverse temperature /3. 
Upper curve: critical point (/? > 3, /? = pc) (see ( 6.25D ). Lower curve: condensed phase (/? > 2, 



p > Pc) (see ( [7.7D ). Thin dashed lines: continuation of the results to high temperature. 



8 Discussion 

At the onset of this paper we gave a comparative presentation of two classes of dynamical 
urn models, the Ehrenfest class and the Monkey class. All these models have simple 
static properties, because their Hamiltonian is a sum of contributions of independent 
boxes. They possess however interesting nonequilibrium dynamical properties, even in 
the mean-field geometry. The backgammon model |, ^, [1^ is a prototypical example 
of the Ehrenfest class, while model B of Reference 0, and the zeta urn model [|ll|, |12[ 
which is the subject of the present work, are examples of the other class. 

Let us come back to the role of statistics in the definition and properties of such 
models. The essential difference between the Ehrenfest class and the Monkey class indeed 
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resides in matters related to a priori statistics. Statistics enters the dynamical definition 
of models: the proposed moves for the Ehrenfest class (respectively, the Monkey class) 
are chosen according to ball-box statistics (respectively, box-box statistics). Consistently, 
statistics also enters their equilibrium definition: the partition functions ( |1.2| ), [ ll-'Sj ] of 
the original Ehrenfest model, and of its generalization to M urns, involve a factorial 
of the total number of particles in their denominators. Inverse factorials l/A^J, taking 
into account equivalent labelings of particles within each box, are also involved in the 
evaluation of the partition function (|1.6|) . 

In statistical mechanics with Maxwell-Boltzmann statistics, the presence of inverse 
factorials has its origin in the indiscernibility of identical classical objects. These facto- 
rials are absent for the Monkey class, in which the populations Ni are involved in flat 
sums, such as ( [1. 71 ), just as occupation numbers are in quantum-mechanical statistical 
mechanics with Bose-Einstein statistics. There is of course nothing quantum-mechanical 
in the urn models considered here. It can be said for short, following Reference [11^1 , that 
the equilibrium statistics of the Ehrenfest class is Maxwell-Boltzmann, while that of the 
Monkey class is Bose-Einstein. Table 1 illustrates this discussion. 



Class 


Ehrenfest 


Monkey 


Dynamical rule 


ball-box 


box-box 


Equilibrium weight of box i 






Statistical-mechanical analogue 


classical 
(Maxwell-Boltzmann) 


quantum- mechanical 
(Bose-Einstein) 



Table 1: Comparison of a priori statistics for the Ehrenfest class and the Monkey class of 
dynamical urn models: dynamical and equilibrium aspects. 



The main focus of this work concerns the nonequilibrium properties of the zeta urn 
model, with emphasis on the aging behavior of the two-time correlation and response 
functions of the fluctuating population of a given box, and of the corresponding fluc- 
tuation-dissipation ratio. We considered successively the critical hne (p = pc) and the 
condensed phase (p > pc). We summarize these two cases below. 

The critical line in the temperature-density plane corresponds to a line of fixed points, 
parametrized by inverse temperature (3 |]TT|. In other words, critical exponents, both 
static and dynamical, and more generally universal quantities, depend continuously on 
temperature. In the present context, a universal quantity ought to be independent: 

- of the initial state (provided it is homogeneous and disordered, i.e., the /a;(0) decay 

rapidly), 

- of the specific form of the energy of each box (provided it diverges logarithmically, as 

E{Ni) InA^j, for large occupation numbers). 
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- and of details of the dynamics (such as Metropohs versus heat-bath). 



Our results ( |6.5| ) and ( |6.20| ) for the two-time correlation and response function have 
the expected product form (see [T^ and references therein), involving: a common non- 
universal prefactor, a negative power of the waiting time, related to the anomalous dimen- 
sion of the observable, and a universal scaling function of the time ratio, or temporal as- 
pect ratio, X = t/s. Accordingly, the fluctuation-dissipation ratio (|6.23|) , X(t, s) ~ X{x), 
only depends on x in the nonequilibrium scaling regime. The limit value Xoo = ^"(00) 
has been recently emphasized fl^, [l^ to be a new universal quantity, characteristic of 
nonequilibrium critical dynamics. For the zeta urn model, the result ( |6.25| ) applies to 
the regular part of the critical line (/? > 3), so that 4/5 < X^o < 1, as shown by the 
upper curve in Figure |I[ This range is unusual for a critical system. Indeed, statistical- 
mechanical models such as ferromagnets are observed to have < < 1/2 at their crit- 
ical point. The upper bound X^o = 1/2, corresponding to the mean-field situation W7 



is also observed in a range of simpler models [13, 14, 0. It is worth noticing that the 
backgammon model also has a high fluctuation-dissipation ratio at its zero-temperature 
critical point, namely X{t, s) ~ 1 — C/(lns)^ for s -C t, where the amplitude C depends 



both on the observable and on the dynamical rule 10|. 



The present analysis of condensation dynamics for p > pc extends and completes 



that begun in |T2[. We have investigated various quantities related to the occupation 
probabilities in the scaling regime, describing the growing condensate. An approximate 
analysis of the Heun operator ( |3.14| ) has allowed us to obtain asymptotic expressions 
for various quantities at high and low temperature. The fluctuation-dissipation ratio 
admits a non-trivial limit value X^o throughout the condensed phase, shown by the lower 
curve in Figure |l|. Expression ( [7. 91) shows that Xqo ~ slowly goes to zero at low 

temperature. This behavior is very different from that of conventional models. It is 
indeed currently accepted pO| that coarsening systems, such as ferromagnets quenched 
from a high-temperature initial state, have identically X^o = throughout their low- 
temperature phase, i.e., for any temperature below T^. These unusual features of the 
zeta urn model are less of a surprise if one remembers that the condensation dynamics 
of the present model is basically different from a domain-growth or coarsening dynamics. 
In the latter case, phase separation takes place in a statistically homogeneous way, at 
least for an infinite system. To the contrary, in the present situation, condensation takes 
place in a very inhomogeneous fashion. The form of equation ( p.9|) indeed demonstrates 
that the condensate of particles is shared by an ever decreasing fraction of boxes, scaling 
as each of them having a population growing as (until size effects eventually 

become important for t ~ JVP). 
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A Solving equation ( |St6i ) by Mellin transformation 

The Mellin transformation provides an efficient alternative way of solving the differential 
equation (|3.(j|). 

We define the Mellin transform Mf{z) of a function f{u) by 

noo f A 2 

Mf{z)= u~^~'f{u)du, f{u)= —u^Mfiz). (A.l) 

JO J 2,711 

Equation ( p.6|) , together with the boundary condition -F(O) = 1, is equivalent to the 
following functional equation 

Mi-f{z + 2) z 



Mi_f{z) 2{z + 2){(3-1-z) 

for the Mellin transform of 1 — F{u). The solution of this functional equation, with 
boundary condition F(oo) = 0, reads 

^^-^^^^ = ^^TFTmY (0<Rez</5 + l). (A.2) 
Similarly, for Re^; < 0, the Mellin transform of F{u), 

^H^) = -^^^Prj (Re.<0), (A.3) 

is the formal opposite of equation (|A.2|) . 

We also give for further reference the expression of the Mellin transform of the deriva- 
tive F'{u): 

MF'iz) = -iz + l)Mi_p(2 + 1) = (z + 1)Mf{z + 1) = ^ {Rez < p). 

^ . ' . ' 2^+ir(^ 

-l<Re2</3 Re2<-1 V ^ / 

(A.4) 



B Solving equation (|5T3|) by separation of variables 
and spectral superposition 

Partial differential equations such as ( |6.3|) , with given initial and boundary conditions, can 
be solved explicitly by the method of separation of variables and spectral superposition. 

To do so, we need a basis of eigenf unctions of the differential operator V of equa- 
tion ( p.7|) . Inspired by the explicit form of the solution ( p.8|) , we set 

G{u) = u^+^ e""'/^ L{v), v = u^A. 
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The eigenvalue equation {V — E)G{u) = becomes 



d^L f(3 + 3 ^ 



dL 

dv 



+ {E-1)L = 0, 



which is the differential equation obeyed by the Laguerre polynomials L"(v) |21|, with 
a = {(3 + l)/2 and n = E - 1. 

The eigenvalues of the operator T> therefore read En = n + 1, with n = 0,1, . . . The 
associated eigenfunctions, 



with 21 



13+3 



-V 



',T(k + ^^) k\{n-k)V 



(B.i; 



obey the orthogonality property 

Gm{u) Gn{u) e"'/^ du = Nn Sm,n, 



(B.2) 



The Mellin transform Mq^ {z) of these eigenfunctions can again be evaluated in closed 
form: 

2/3-r(n + ^) r(^ 



Mniz) 



(Rez <f3 + l), 



(B.3) 



where the normalization has been fixed by the condition L^^+^)/^(0) = 1. 

The method of spectral superposition consists in looking for the solution G{u,x) of 
equation ( |6.3| ) as a linear superposition of the form 

oo 



u) 



n=0 



The coefficients a„ are determined by the initial condition 
property ( |B.2D , we obtain 

an = ^ ru-F{u)Lr'y\u^/A)du. 

Nn Jo 

Then, using (p.l| ) and (|B.3|) , we are left with 



(B.4) 

Using the orthogonality 



9I-/3 n 



^^{2k + 3)k\{n-k)\T[k + 



§+3 
2 



(B.5) 
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